Mode-coupling theory for structural and conformational dynamics of polymer melts 
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A mode-coupling theory for dense polymeric systems is developed which unifyingly incorporates 
the segmental cage effect relevant for structural slowing down and polymer chain conformational 
degrees of freedom. An ideal glass transition of polymer melts is predicted which becomes molecular- 
weight independent for large molecules. The theory provides a microscopic justification for the use 
of the Rouse theory in polymer melts, and the results for Rouse-mode correlators and mean-squared 
displacements are in good agreement with computer simulation results. 
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Polymeric materials often find applications as amor- 
phous solids or glasses. During their processing, the 
structural relaxation in polymeric melts plays an im- 
portant role because it affects transport and viscoelastic 
properties. An understanding of conformational dynam- 
ics is also required for single-chain diffusional processes. 
Whereas a number of phenomenological concepts have 
been invented, no microscopic theory exists which suc- 
cessfully explains the slow structural and conformational 
relaxations of polymeric melts close to the glass transi- 
tion, as for example studied in detailed computer simu- 
lations of model systems [Q . 

We present a first-principles theory for structural and 
conformational dynamics of uncntangled polymer chains 
using the site formalism which has already provided a 
good description of their static packing Building 
upon the mode-coupling theory (MCT) for simple and 
low-molecular systems ^ |^, our macromolecular ex- 
tension captures the segmental "cage effect" which causes 
a (asymptotically) molecular-weight-independent ideal 
glass transition driven by steric hindrance on microscopic 
length scales. A decoupling of local collective structural 
relaxation from single-chain conformational fluctuations 
is obtained, leading to a first-principles derivation of the 
Rouse model Q . This widely used model assumes a chain 
to be in a Markovian heat bath, and therefore its appli- 
cability requires a time-scale decoupling as our theory 
achieves. Entanglement effects are neglected in our 
work because we do not ensure uncrossability of chains 
and consider only local isotropic forces 0. 

A system of n chains, consisting of N identical 
monomers or segments, distributed with density p is con- 
sidered. In the site formalism the structural variables 
are the monomer-density fluctuations for wave vector q, 
p't = ^"^j^ exp(i(f • r.°-), where r° denotes the position 
of ath monomer in ith chain. Structural dynamics shall 
be described by coherent density correlators, F^^{t) = 
{p°j[t)* P^^ In. Here (• • • ) denotes canonical averaging for 
temperature T. A single (or tagged) polymer (labeled s) 
exhibits density fluctuations p°j ^ = exp(zg • r°), and cor- 
relators, F^^^(t) = (p|s(t)*p^s), characterize the single- 



chain dynamics. 

We apply MCT equations in the site representation |^ 
to flexible macromolecules expecting them to capture in- 
termolecular caging in polymeric melts. The required 
inputs are the static structure factors S^^ = F°^(0), 
^ab _ pab^ j-q-j ^ g^j^^ ^-^^ dircct Correlation functions c"'' [D . 
There are severe difficulties to solve the N x A^-matrix 
MCT equations for polymers because the degrees of poly- 
merization N of interest are large. The simplifications 
adopted here are to neglect chain-end effects for 

= Cq, and to consider the site-averaged correlator, 

Fq{t) = (1/iV) Y.a,b=i Fq^t), which deals with the total 
monomer-density fluctuations. This mean-field like ap- 
proximation replaces the site-specific surroundings of a 
segment by an averaged one. It is supported by the ob- 
servation that Sq = Fq{0) captures the static correlations 
on the segmental length scale |^ . 

We find a set of scalar equations for the normalized 
coherent correlator cj)q{t) — Fq(t)ISq: 

dU,{t) + nl<Pq{t) 

+ [ dt'mqit^t')dt'(j}qit')=0, (1) 

Jo 

mq{t) dkV{q; k,p) Mt) Mt)- (2) 

Here fl^ = q^v^lSq and V = PmSqSkSp{q ■ [kck + 
pcp]}^/(27r)'^g'* with v denoting the monomer thermal ve- 
locity, pm = Np, and p = q — k. These equations are 
formally identical to the ones for simple systems Q. On 
the other hand, one finds for the single-chain dynamics: 

d!Fq.s{t)+nl,FqAt) 

+ nlf dt'inq,s{t~t')dt'Fq,S')^0, (3) 

JQ 

yjac I- 

<s{t) = E ^ / '^^^^('?*; ^'P*) w Mt)- (4) 

Here ^ = ^2^,2^-1 ^nd = Pm{q-plqf Spcll{2TTf. 

A traditional description of the single-chain dynam- 
ics is in terms of Rouse modes [0. For discrete chains. 
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FIG. 1: The solid line exhibits the glass-form factors at 
ifi = ipc for A'' = 10 along with the simulation result (circles) 
for a similar system taken from Ref. |l4[ The dashed line 
denotes Sq multiplied by 0.15. The inset shows ifidN) rescaled 
by vjj;"' = ipciN ^ oo) « 1.309 as a function of log^Q N. 



Rouse-mode correlators are defined by Cp,p'{t) — {Xp{t) ■ 
Xp')/3N with Xp = v^Ef=i ^7 cos[(a - 1/2)ptt/N]. 
Since ~ l+iq-f^ for small q, Cpp' (t) can be expressed 
as a linear combination of F^\{t) for g ^ 0. Here, it is 
neither assumed that Cppi (t) are diagonal, nor that the 
decay is exponential. Therefore, our microscopic theory 
can test the validity of the Rouse theory. 

The polymers considered here shall be modeled as fol- 
lows. First, chains are assumed to be Gaussian, for 
which w^'^ — exp[— q^|a — c|&^/6] with the statistical 
segment length b. Second, each monomer is modeled 
as a hard sphere of diameter d and mass m, and we 
set b — d. All equilibrium properties are then speci- 
fied by the packing fraction ip — Trp^d^ /Q and the degree 
of polymerization N. Third, Sq and Cq are evaluated 
from the polymer RISM theory Let us add that the 
value of if can become larger than 1 for very dense sys- 
tems since corrections due to unphysical intrapolymer 
monomer overlap |^ are not taken into account in the 
present work. From here on, the units will be chosen so 
that d = V = m = 1. 

The derived MCT equations (|l|) and (||) exhibit an 
ideal liquid-glass transition upon increasing (/? at a criti- 
cal value ipc as discussed for simple systems |l^ . The 
solid line in Fig. |l| exhibits the resulting glass-form factors 
fq ~ 't'qi't' ^ co) at Lpc for N = 10, and the dashed line 
Sq multiplied by 0.15. The fq also measure the ampli- 
tude of the slow structural relaxation in the liquid state. 
Inspection of Eq. (||) shows that the liquid-glass transi- 



tion is driven mainly by the changes of Sq for gmax ~ 7 
where is maximal, i.e., by local fluctuations connected 
to the average nearest-neighbor-monomer distance (the 
cage effect). With varying g, oscillates in phase with 
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FIG. 2: Rouse-mode correlators Cp{t) of p = 1 (solid), p — 5 
(dashed), and p = 9 (dotted lines) for A'' = 10 and {ip — 
'Pc)/'Pc ~ —10"'^. The horizontal lines mark the plateaus 
fp. The inset shows the ratio Tp/rqmax of the a-relaxation 
times (circles) and the stretching exponent /3p (squares) as a 
function of p. The dotted line denotes a fitting based on the 
formula Tp = C/12 sin^(p7r/27V). 



for a similar system 0, |T|, |T|, |T|l. 

semiquantitative especially for q « 



The agreement is 
?mav, which is the 



relevant q range for the ideal glass transition ||15|. The 
inset demonstrates that ipc{N) becomes independent of 
N for large N. This is because the glassy arrest is driven 
by the local cage effect, and the global chain size plays a 
gradually smaller role for larger N . 

For packing fractions tp close to but below pc, corre- 
lator (j)Ait) of any variable A coupling to density fluc- 
tuations exhibits slow structural relaxation and decays 
in two steps: the decay towards the plateau fi, fol- 
lowed by the decay to zero (so-called a-decay) Q. 
Detailed analysis of (l)q{t) and F^\{t) will be presented 
in later publications. Here, only the Rouse-mode cor- 
relators Cpp'{t) shall be studied. Typically, their off- 
diagonal elements are found to be much smaller (at 
most only a few %) compared to the diagonal ones at 
all times. Therefore, the normalized diagonal elements, 
Cp{t) — Cpp{t)/Cpp{0), shall be considered. Figure || ex- 
hibits representative results for iV = 10 for a reduced 
packing fraction (ip — ipc)/p^c = —10^^. 

The curves shown in Fig. || do not clearly exhibit 
the two-step-relaxation scenario. This is because the 
plateaus of Cp{t) are so large, fp > 0.9, that only less 
than 10% of the decay are left for the relaxation towards 
the plateau. Thus, most of the relaxation of Cp(t) occurs 
in the a regime. Let us characterize the a-time scale by 
Cp(Tp) = fp/'^O- The corresponding a-time scale Tqmax 
shall be introduced for the coherent correlator (f>q{t) for 
q = Qn 



characterizing the local dynamics of the sur- 
Sq as found in simple systems |l^. The circles denote rounding medium. The ratio Tp/rqmax is shown in the 
the result of molecular-dynamics simulation performed inset as a function of p. For small p, the scales Tp are 
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much larger than the local scale Tqmax- 

The Rouse theory assumes that all dynamical correla- 
tions in the surroundings are much faster than the single- 
chain dynamics |^]. Since a polymer is surrounded by 
identical polymers, the assumption of the time-scale sep- 
aration cannot be justified a priori. Our microscopic the- 
ory verifies this central assumption. The Rouse theory 
predicts, within our units, Tp = C/12 sin^(p7r/2iV) where 

denotes the monomer friction. This formula, using C 
as a fit parameter, is shown as the dotted line, and we 

found C/T-qmax = 16.5. 

The shape of Cp{t) in the a regime is often character- 
ized by the stretching exponent f3p of the Kohlrausch-law 
fit: Cp{t) oc exp[-(t/T;)'''']. We found /3p > 0.9 for all p 
as shown in the inset. For small p, (ip is close to 1 due 
to the large separation of the scales Tp and Tqmax as dis- 
cussed above. /3p decreases as p increases since for large p 
the scales Tp become comparable to Tqmax- However, I3p 
remains close to 1 even for large p. This is because their 
plateaus are high as explained in Ref. Thus, all 
our Cp{t) exhibit nearly Debye-relaxation as assumed by 
the Rouse theory. Let us note that the found features for 
Cp{t) hold also for N = 100, and are in agreement with 
simulation results . 

Of interest are the monomer mean-squared displace- 
ment (MSD) ga{t) = ([rf(0-rf(0)]2) and the center-of- 
mass MSD, gc{t)- Let us also introduce the monomer- 
averaged MSD, cjMit) = {l/N)j:^^^ga{t). Figure I 
exhibits gM{t) and gc{t) for TV = 10 and 100 for 
(ip — ipc)/^c = — 10~^. For short times, the MSD ex- 
hibits ballistic motion: gx{t) (x {X = M or C). As the 
time increases, the MSD begins to be suppressed due to 
the cage effect, and there appears the so-called /3 regime 
where gx{t) is close to the square of the critical local- 
ization length, gx{t) « 6(r3f)^. For this regime, there 
holds I, I |l| 

gx{t) = 6(r^)2 - 6hxG{t), X - M, C, (5) 

where G{t) denotes the (3 correlator and hx the critical 
amplitude. Thus, gmit) and gc{t) cross their plateaus 
6(r^)^ at the same time. The value r^j « 0.077, quan- 
tifying the monomer localization, for both iV = 10 and 
100 is consistent with Lindemann's melting criterion, 
is reduced by about \/\/N compared to r^^, as expected 
for independent motions of constituent monomers. The 
asymptotic law (|^) for each MSD is drawn as dashed 
line, and its range of validity is indicated by diamonds. 
The trends in gx (t) up to the /3 regime are qualitatively 
similar to the ones for a sphere in a simple system . 

The increase of gx{t) above the plateau 6(r^)^ to- 
wards the diffusion asymptote, gx(t) ~ 6Dt with the 
diffusivity D, is the a process of the MSD. In contrast 
to gc{t), gM{t) in this regime is significantly affected by 
chain connectivity since the segments participate in the 
conformational motion and most of the relaxation of Cp(t) 




-2 2 4 6 8 




FIG. 3; Double logarithmic presentation of the MSDs gM{t) 
and gcit) (solid lines) for A'' = 10 (upper) and 100 (lower 
panel) for {(p — (pc)/ipc ~ — 10~^. The dashed lines denote 
the asymptotic formula given in Eq. (^. The open diamonds 
mark the points where the dashed lines differ from the solid 
ones by 10%. The dashed-dotted lines exhibit the diffusion 
law, 6Dt. The dotted lines show the power laws, ~ t^, 
with X specified in each panel. The horizontal dotted lines 
from bottom to top successively mark the plateau 6{rc)^ for 
gc{t), the plateau 6(r^)^ for gM{t), the radius of gyration 
i?3, and the end-to-end distance B? . The inset exhibits the 
ratio gi{t)/gf^/2it)- 



occurs here as explained above. As a result, there ap- 
pears a subdiffusive (^ t^) regime in gM{t)- For N — 10, 
the exponent x is 0.60, which is close to the value (0.63) 
found in the mentioned simulations ||l^. For N — 100, 
we find x = 0.5 as predicted by the Rouse model and by 
the asymptotic evaluation of our theory, and the different 
value for = 10 can be attributed to finite A^ effects. 
Thus, we find a strong polymer-specific effect for the be- 
ginning of the a process, while no such effect is reflected 
in gAiit) up to the P process. 

The insets of Fig. ^ show the ratio of the MSD for end 
and central monomers, gi{t) / gN/2{t)- The ratio starts 
from 1 in the ballistic regime, exhibits a maximum for in- 
termediate times, and tends to 1 in the diffusion regime. 
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It is seen that the ratio remains close to 1 also within the 
(3 regime. This is because the dynamics here does not 
reflect the chain connectivity as explained above. The 
Rouse theory predicts the maximum to be 2 within the 
time regime where the monomer MSD exhibits the t^^"^ 
law, i.e. in the a regime. The result for N = 100 clearly 
indicates this behavior, while the maximum for A'^ = 10 
is slightly smaller (1.8) due to finite N effects. The shape 
of the ratio for = 10 in the (3 and a regimes is in semi- 
quantitative agreement with the mentioned simulation 
results [19 . Thus, the chain-end effect is properly taken 



into account in our theory even with neglecting that effect 
for Cq^. This is because the matrix structure of Eqs. (|^) 
and are preserved for describing the single-polymer 
dynamics. 

In summary, a first-principles theory for structural 
slowing down of dense polymer systems has been pre- 
sented which also provides a microscopic derivation of 
the Rouse model for unentangled chain melts. Chain 
connectivity is seen to cause polymer specific long-time 
anomalies of the a-process. We use the concept of an 
ideal MCT glass transition, familiar for simple liquids 
and colloidal suspensions and our results agree well 
with simulation studies of simple (coarse grained) poly- 
mer models. Our theory thus provides insights into dy- 
namical aspects typical for polymer melts which are also 
observed in models j2^, ^ with a more realistic local 
chemistry. 

We thank J. Baschnagel, M. Aichele, and W. Gotze 
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